Performance comparison of dynamical decoupling sequences for a qubit in a rapidly 

fluctuating spin-bath 



Gonzalo A. Alvarez, 1 >0 Ashok Ajoy, 1,2,3 Xinhua Peng, 1,4 and Dieter Suter^Q 

1 Fakultat Physik, Technische Universitat Dortmund, D-44221 Dortmund, Germany. 
2 Birla Institute of Technology and Science - Pilani, Zuarinagar, Goa - 403726, India. 
3 NMR Research Centre, Indian Institute of Science, Bangalore - 560012, India. 
Hefei National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, 
University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China 

Avoiding the loss of coherence of quantum mechanical states is an important prerequisite for quan- 
tum information processing. Dynamical decoupling (DD) is one of the most effective experimental 
methods for maintaining coherence, especially when one can access only the qubit-system and not 
its environment (bath). It involves the application of pulses to the system whose net effect is a 
reversal of the system-environment interaction. In any real system, however, the environment is 
not static, and therefore the reversal of the system-environment interaction becomes imperfect if 
the spacing between refocusing pulses becomes comparable to or longer than the correlation time 
of the environment. The efficiency of the refocusing improves therefore if the spacing between the 
pulses is reduced. Here, we quantify the efficiency of different DD sequences in preserving different 
quantum states. We use 13 C nuclear spins as qubits and an environment of 1 H nuclear spins as the 
environment, which couples to the qubit via magnetic dipole-dipole couplings. Strong dipole-dipole 
couplings between the proton spins result in a rapidly fluctuating environment with a correlation 
time of the order of 100 /is. Our experimental results show that short delays between the pulses 
yield better performance if they are compared with the bath correlation time. However, as the pulse 
spacing becomes shorter than the bath correlation time, an optimum is reached. For even shorter 
delays, the pulse imperfections dominate over the decoherence losses and cause the quantum state 
to decay. 

PACS numbers: 03.65.Yz,03.67.Pp,76.60.-k ,76.60.Lz 



I. INTRODUCTION 

Quantum mechanical systems have an enormous po- 
tential for realizing information processing devices that 
are qualitatively more powerful than systems based on 
classical physics The main requirement for quantum 
information processing (QIP) is that the system evolves 
according to the Schrodinger equation, under the influ- 
ence of a Hamiltonian that is under precise experimental 
control. However, no system is completely isolated, and 
disturbances from its surrounding environment (bath) 
spoil the quantum identity of the system. This process is 
often called decoherence [2] and limits the time scale over 
which quantum information can be retained and the dis- 
tance over which it can be transmitted Reducing 
the effects of decoherence is therefore one of the main re- 
quirements for reliable quantum information processing. 
Several protocols have been developed for quantum error 
correction 0, 0|; however, they prove advantageous only 
for low levels of environmental noise. 

In all existing experimental architectures for QIP the 
noise background is too large, and this limits the appli- 
cability of these protocols. A promising technique for 
reducing the noise to a level where error-correcting codes 
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can take over is called Dynamical Decoupling (DD) 
It aims to reduce the interaction of the system with the 
environment through control operations acting only on 
the system. It requires relatively modest resources, since 
it requires no overhead of information encoding, measure- 
ments or feedback. 

Although the mathematical framework of dynamical 
decoupling was introduced fairly recently [81, t he six- 
decade old Hahn NMR spin-echo experiment can be 
considered as the earliest and simplest implementation of 
this method. It consists of the application of a 7r-pulse to 
a spin qubit ensemble, at time r after the spins were left 
to undergo Larmor precession in a magnetic field. This 
effectively reverses a pure dephasing system-environment 
(SE) interaction, i.e. one that does not cause a net ex- 
change of energy between the system and the bath. The 
combined effect of the evolution before the refocusing 
pulse and a second period of the same duration after the 
pulse vanishes. Physically, the dephasing and rephasing 
of the spins can be observed as an apparent decay of the 
average magnetization in the system and a subsequent 
increase after the refocusing pulse (a spin echo). 

The Hahn-echo can (for an ideal pulse) completely 
eliminate the interaction with the environment, provided 
it is time-invariant. In practice, this is often not the case, 
and a change in the environment reduces the refocusing 
efficiency flix [Tl| . To reduce the problems due to a time- 
dependent environment, Carr and Purcell suggested to 
replace the single pulse of the Hahn echo by a sequence 
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of pulses at shorter intervals (the CP sequence) , thus 
reducing the changes in the environment between succes- 
sive pulses. For sufficiently short pulse intervals, elimina- 
tion of system-environment interactions became possible 
even in a time-dependent environment. However, the in- 
creased number of pulses led to another problem: if the 
refocusing pulses are not perfect, they actually become 
a source of decoherence (and thus signal loss) instead of 
eliminating it. This problem was significantly reduced 
by a simple modification of the CP sequence: if the rota- 
tion axis of the refocusing pulses is parallel to the initial 
spin orientation, the effect of pulse errors is significantly 
reduced over a cycle [l2l |. This is known in literature as 
the CPMG sequence. 

In the context of QIP, there has been renewed effort in 
eliminating the effects of the system-environment inter- 
action that lead to the loss of quantum information. For 
this, it is often important that the refocusing reduces the 
effect of the system-environment interactions by several 
orders of magnitude. In addition, the effect of pulse er- 
rors must be minimized, and the sequence has to work 
for all possible initial states of the system. Several pulse 
sequences that achieve this were introduced fl3l - [l5j |. 
which consist of periodic sequences of pulses; they were 
thus called periodic DD (PDD). By design, they allow 
one to decouple the system from the environment for a 
general SE interaction, i.e. one that causes dephasing as 
well as dissipation. 

Experimentally, DD is achieved by iteratively applying 
to the system a series of stroboscopic control pulses in 
cycles of period r c . Over that period, the time-averaged 
SE interaction Hamiltonian vanishes. The time average 
over T r can be calculated using average Hamiltonian the- 
ory [l6| ■ If the average Hamiltonians are calculated by a 
series expansion, such as the Magnus expansion, improv- 
ing the pulse sequence usually corresponds to progres- 
sively eliminating higher order terms in the expansion. 
Khodjasteh and Lidar [13] introduced concatenated DD 
(CDD) as a scheme that recursively generates higher or- 
der DD sequences for this purpose. Here, the lowest level 
of concatenation is a PDD sequence. The improvement 
achieved by concatenation comes at the expense of an 
exponential growth [l7| in the number of applied control 
pulses. In contrast, for the case of a pure dephasing or 
pure dissipative interaction Hamiltonian fl8L fl9l | , Uhrig 
developed a sequence (UDD) [20] that reduces higher or- 
ders in the Magnus expansion with only a linear overhead 
in the number of pulses. Unlike other DD sequences, in 
the UDD sequence the delay between successive pulses 
is not equal, i.e. the pulses are not equidistant. In the 
limit of a two-pulse cycle, UDD reduces to the CPMG 
sequence. Recent proposals of DD sequences that are a 
hybrid between UDD and CDD are predicted to improve 
DD performance of previous methods [2ll |2^. 

The UDD sequence was tested on ion traps (23l . [2~3 | , 
electron paramagnetic resonance [25| and liquid-state 
NMR [26], and found to outperform equidistant pulse 
sequences, in particular CPMG, for environments with 



a high-frequency or strong cutoff. CDD sequences were 
recently tested in solid-state NMR [27} ■ However, while 
some sequences for particular environmental noises were 
tested, a comparison between sequences for different 
kinds of environments is still missing. Most of the se- 
quences were designed assuming ideal pulses and some 
of them predict to compensate pulse imperfection. How- 
ever an experimental test of this aspect is still needed. 
In parallel to this work recent DD implementations on 
a qubit interac ting with a slowly fluctuating spin-bath 
were tested [28M32| . 

Other questions relate to the optimal cycle time: It 
is theoretically predicted and experimentally demon- 
strated that sequences that reduce higher order terms 
of the Magnus expansion perform better than low or- 
der sequences for slow motion environments with high- 
frequency or strong cutoff, when the bath correlation 
time tb is longer than the sequence time r c (cycle time) . 
However the strength and duration of control-pulses are 
limited by hardware, yielding a minimum for the achiev- 
able DD cycle time. As some examples on this direction, 
Viola and Knill proposed a general method for DD with 
bounded controls [15j . Khodjasteh and Lidar, keeping 
the delay between pulses constant, predicted an optimal 
CDD order for reducing decoherence |33l. l34j. Biercuk et 
al. [23l l2~j | needed to consider the finite length of pulses 
in their simulations, assuming them perfect but produc- 
ing a spin-lock during their application times, in order to 
fit them to the experiments. Additionally Hodgson et al. 
[35[, while assuming instantaneous perfect pulses, theo- 
retically analyzed DD performance constraining also the 
delay between pulses. They set a lower limit to the delays 
making them larger than the pulse duration to satisfy 
the instantaneous pulse approximation in their theoreti- 
cal model for experimental conditions. When the delays 
are strongly constrained, they predict that DD protocols 
like CDD or UDD, which are designed to improve the per- 
formance of lower DD orders if the regime of arbitrarily 
small pulse separations is achievable, in general lose their 
advantages. However, experiments are missing in order 
to demonstrate these predictions and very little is known 
about the performance of DD sequences under conditions 
where the cycle times are comparable to or longer than 
the bath correlation times. Recently Pryadko and Quiroz 
approached this regime, but only for the extreme case of 
a Markovian environment [36| . 

While the finite length of pulses limits the minimum 
cycle time reducing the maximal achievable DD perfor- 
mance, their imperfections also contribute to reducing 
it. It is well know that CPMG-like sequences are too 
sensitive to the initial state when pulse errors are con- 
sidered fl3l . [l4| . A comparisons between the CPMG and 
UDD sensitivity against pulse errors was performed in 
Ref. (2^. Overall UDD was shown to be more robust 
against flip angle errors and static offset errors, with the 
exception that CPMG is more robust for initial states 
longitudinal to the control pulses. But both of them are 
too asymmetric against initial state directions. In gen- 
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eral, while some DD sequences were developed to com- 
pensate flip-angle errors and to have a performance more 
symmetric against initial conditions, an extensive study 
of their performance from a QIP perspective is still miss- 
ing and additionally is not done for CDD sequences. For 
example, an optimal cycle time when considering imper- 
fect finite pulses was predicted by Khodjasteh and Lidar 

In this article, we compare experimentally the perfor- 
mance of different DD sequences on a spin-based solid- 
state system where the cycle time t c is comparable to or 
longer than the correlation time tb of the environment. 
Here, the spin(qubit)-system interacts with a spin-bath 
where the spectral density of the bath is given by a nor- 
mal (Gaussian) distribution. This kind of systems, typi- 
cal in NMR [37j , are encountered in a wide range of solid- 
stat e sy stems, as for example electron spins in diamonds 
[30l - [32l |. electron spins in quantum dots [HI, [2j| |38| and 
donors in silicon (39l |40| which appear to be promising 
candidates for future QIP implementations. In particu- 
lar we consider the case where the interaction with the 
bath is weak compared with the intra-bath interaction. 
For one side the latter point complement and distinguish 
our work from the recent submitted articles (28l - l32l |. For 
the other the aim of our work is a comprehensively and 
detailed comparison of the performance of different se- 
quences considering different initial states. We find how 
the performance of the DD sequences depends on the 
initial state of the qubit ensemble with respect to the ro- 
tation axis of control pulses with finite precision. When 
they are in the same direction, the CPMG sequence is 
the best DD sequence for reducing decoherence - i.e. 
it maintains the state of the ensemble for the longest 
time. However, if the initial state of the ensemble is not 
known, we find that concatenated dynamical decoupling 
(CDD) provides the best overall performance. Stated 
equivalently, the CDD scheme provides the best over- 
all minimization of the environmentally driven quantum 
mechanical evolution of the system. Additionally we ex- 
perimentally demonstrate and quantify the predicted op- 
timal delay times for maximizing the performance of the 
respective DD sequences. This implies that pulse errors 
are a limiting factor that must be reduced to improve DD 
performances. In general our results complement some 
of the previous findings and predictions for some of the 
experimentally tested DD sequences and provide new re- 
sults for untested ones. One of the main message is that 
a fair comparison of the performance of DD sequences 
should use a constant average number of pulses per unit 
time. 

This paper is organized as follows. Section ITI1 describes 
the qubit and bath system used in our experiment, and 
the mechanisms of coupling between them. In section ITTT1 
we give a brief summary of dynamical decoupling and 
a description of the tested sequences - the Halm Echo, 
CPMG, PDD, CDD and UDD. Our limited choice of se- 
quences includes those most accepted by the QIP commu- 
nity and allows us to discuss the most important points. 



Section IIVI contains the experimental results and their 
analysis. In section [V] we compare the various DD se- 
quences under the same conditions. In the last section 
we draw some conclusions. 



II. THE SYSTEM 

Our system consists of a spin 1/2 (qubit) in a strong 
magnetic field oriented along the z-axis, interacting with 
a bath consisting of a different type of spins 1/2. The 
total Hamiltonian in the laboratory frame is 



n L = U h a 



U 



SE 



Ei 



(1) 



where Jig is the system Hamiltonian, "H^ is the environ- 
ment Hamiltonian and Hg E is the system-environment 
interaction Hamiltonian: 



Us = u S S z 



U L SE = S z J2b S] H, 
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(2) 
(3) 

21*%- {rJi + PJi)] , (4) 



i<3 



where S is the spin operator of the system qubit, the spin 
operators i|, P y and J| act on the j th bath spin, ojs and 
loi are the Zeeman frequencies of the system spin and 
the bath spins respectively, bsj and dij are the coupling 
constants, and we use frequency units (H = 1). In solids, 
the spin-spin interaction is dominated by the dipolar in- 
teraction [37| . Since S and / are different types of nuclei, 
it is possible to neglect the terms of the dipolar coupling 
Hamiltonian that do not commute with the strong Zee- 
man interaction because \bsj\ / \u>s — "C 10~ 4 [37| . 
The remaining terms have the Ising form Similarly, 
the homonuclear interaction between the bath spins is 
truncated to those terms that commute with the total 
Zeeman coupling, which we assume to be identical for all 
bath spins. 

In the high-temperature thermal equilibrium (37^ . the 
density operator of the system spin is 



PS,eq. oc S z 



(3) 



where we consider only the system (S'-spin) part of the 
total Hilbert space. We also neglect the part proportional 
to the unit operator, which does not evolve in time and 
does not contribute to the observable signal. 

To generate the initial state for our DD measurements, 
we rotate the thermal state to the xy-plane by applying 
a 7r/2 pulse. The resulting state is 



ps(0) ex S 



{;}■ 



(6) 



For an isolated spin system this magnetization precesses 
indefinitely around the static magnetic field at the Zee- 
man frequency u>s- 
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Taking the system-environment interaction into ac- 
count, the effect of the coupling operator H$ E is the gen- 
eration of product terms of the form S±P Z in the density 
operator, correlating the system with the environment. 
Since we only observe the system part of the total Hilbert 
space, we effectively project the correlated system onto 
this subspace, 



Ps = Tr/ {ptot} ■ 



(7) 



where Tr/ represents the partial trace over the environ- 
mental degrees of freedom and ptot represents the den- 
sity operator of system plus environment. The result 
of this projection corresponds to a loss of coherence by 
dephasing. This free evolution of the system under the 
SE interaction is called the free induction decay (FID) 
in NMR terminology. The decay process of the system 
state is usually called relaxation in NMR terminology, or 
decoherence in quantum information. 

In the following, we will describe the dynamics of the 
system in a rotating frame of reference [33|: The system 
rotates at the (angular) frequency ajg around the z-axis 
and the environment at wj. As a result, the rotating 
frame Hamiltonian becomes 



where 



H f = H s + H S e+He, 



Hs — Hg — uigS z — 0, 
Use = S z ^2b S jI J z , 



(8) 

(9) 
(10) 



Hi 
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This transformation is exact, since the Zeeman terms 
commute with all other terms in the Hamiltonian as well 
as with the equilibrium density operator. 

The effect of the environment-Hamiltonian "He on the 
evolution of the system may be discussed in an interac- 
tion representation with respect to the evolution of the 
isolated environment: the system-environment Hamilto- 
nian then becomes 



fig) (t) = e-^HsE^ 



S,e 



-m E t 



\bs 3 P z 



jn E t 



(12) 



Since He does not commute with Use, the effec- 
tive system-environment interaction Hgg becomes time- 
dependent: the system experiences a coupling to the 
environment that fluctuates. The correlation time 



tb of the time-dependent spin-bath operators P z (t) — 
e -iU E tj^ e iU E t j g g! ennec i by the decay to 1/e of the cor- 
relation function 



Tr 



Tr 



{p z (o)H(t)} 

{imim} 



(13) 



Considering that all the bath spins are equivalent, the 
latter correlation funtions are identical for every j and 
share the same correlation time tb- 



III. DYNAMICAL DECOUPLING 

A. Notation 

The aim of dynamical decoupling is the reduction of 
the interaction of the qubit system with the environment, 
thus retaining the quantum information for as long as 
possible. In the context of DD, it is assumed that it 
is possible to apply arbitrary single-qubit operations to 
the system qubit, but that it is not possible to control 
the environment. One thus applies to the system short, 
strong pulses, whose effect can be described as a refocus- 
ing of the system-environment interaction by the control 

Hamiltonians Hq(S) {t) BUI- 

Let us refer to Fig. Q] and consider a single cycle of the 
sequence having a period r c . In the rotating frame, the 
operator that describes the evolution of the total system 
from to t c is 

N 

U (r c ) = Uf (tat+x) J] Ub (r p ) Of {n) , (14) 

i=l 

where from Eq. © the free evolution operator is 

U f (t) = exp{-iH f t} (15) 

and the control evolution operators that act during the 
time t„ is 



Texp 



dt 



'(H f 



H 



C(S)l 



with T the Dyson time-ordering operator [4l|, |42j. We 
assume that the free evolution Hamiltonian is constant, 
while the control Hamiltonian ^c(S) (*) ^ s cons tant during 
t p but changes for different i. The delay times between 
the control Hamiltonians are Tj = ti — (U—i + t p ) for i = 
2, ..,N + 1 and r\ = t\ — to, where to — 0, tjv+i = t c , and 
U represents the time at which the i th control operation 
starts. Figure [1] shows a graphical representation of these 
definitions. 

Like any unitary evolution, the total propagator can 
be written as the exponential of a Hermitian operator, 



U{t) 



-m e 



(17) 
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Figure 1: Schematic representation of dynamical decoupling 
The solid boxes represents the control pulses . 



Using average Hamiltonian theory [l6| we can calculate 
the effective Hamiltonian Heff as a series expansion, 



The total evolution operator is thus 

N 

U (t c ) = U' fN+1 (t n+1 ,t p ) J] W c (0) U' fi (Ti,r p ) , (22) 

i=l 

where the evolution operators 

U' ft (n,r p ) = Uf ( n ) exp {-iS ei 6i,e} (23) 
represent a modified free evolution. Note that 

U' h [T i ,T p ) = U f {T l ). 

The zero order average Hamiltonian of the free evolu- 
tion periods (|23[) for non-perfect pulses is equivalent to 
interactions of the general form 



Heff = ti [Q) + n {1) + n 



(2) 



(») 



(18) 



The zero order term w-^ is given by the time integral 
of the total Hamiltonian from time to r c . An ideal 
DD sequence makes = %e , i.e. for ideal pulses, 

the interaction Hamiltonian vanishes to zeroth order. In 
the Magnus expansion [43], higher order terms are pro- 
portional to increasing powers of t c /tb, since we assume 
that the environment is weakly coupled to the system 
(pSj TB <C 1) and in consequence tb is the dominant time- 
scale [33] . 

If the basic cycle is iterated M times (see Fig. Q]) , the 
total evolution operator becomes 



U(t = Mr c )= U(t c ) 



(19) 



^■SE — a xS x + a ySy + 0, Z S Z 



= E S^^u + b^H), (24) 

u—x,y,z j 

where a u and Vgj give the renormalized offsets and cou- 
plings respectively, which include the errors of the con- 
trol Hamiltonians. This picture can also consider errors 
of control pulses when r p is comparable with the inverse 
couplings of the free evolution Hamiltonian. 

We now discuss some DD schemes that refocus the 
system-environment interaction. In all these cases, we 
assume that the system is initially prepared in a coherent 
superposition of the computational basis states. We will 
refer to the initial state of the qubit as S x , S y or S z , as 
shown in Fig. [2][a). 



B. Ideal and real pulses 

The usual approximation of hard pulses - having a 
radio-frequency field uj p 3> bsj and duration r p <C 
d^j , bgj implies that we can neglect the free precession 
Hamiltonian and Eq. (fl"6|) simplifies to 



U, 



c (t p ) = exp j-L9, 



(20) 



in the rotating frame, where u = x,y, z and 8 p = lo p t p is 
the rotation angle around the u axis. In what follows, we 
shall denote perfect instantaneous 7r-pulses along x and 

y by X = exp|— iiS x 7r| and Y = exp |— iS y ir^ respec- 
tively, and a free evolution of duration r by f T . 

To take the effect of non-ideal pulses into account, one 
needs to consider errors in the axis and angle of rotation. 
We write the resulting control propagator as the prod- 
uct of the ideal pulse rotation times an error rotation 

<{-iS ei 6 ite \: 
Uh (t p ) = exp {-iS e Ae} expj-i&^p} . (21) 



exp • 



C. Hahn echo 

The Hahn spin-echo experiment (lp| is the pioneer dy- 
namical decoupling method and the building block for 
newer DD proposals. It consists of the application of a 
7r-pulse to the S spin along an axis (say y) transverse to 
the static field Bo at time t causing an echo at time 2t 
[Fig. (Hb)]. The total evolution operator can be summa- 
rized as f T Yf T where the total time (assuming a delta- 
function pulse) is It. As a consequence, the zero-order 
average Hamiltonian is, 



n 



(o) 

Hahr. 



2t 



— j dt'H (t') 



tH S e - tT-Lse 



2t 



0. 

(25) 

The resulting system evolution operator approaches the 

identity to within O (^(t c /tb) 2 ^ . Thus if t c <C Tb, a 

perfect echo (time reversion) is achieved at the total evo- 
lution time t = 2t = t c . When r c is comparable to or 
longer than Tb, the echo decays due to the higher order 
terms. 
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Figure 2: Schemes of dynamical decoupling pulse sequences. 
Empty and solid rectangles represent 7r/2 and 7r pulses re- 
spectively. M represents the number of iterations of the cy- 
cle, (a) Initial state preparation before application of the DD 
sequence, (b) Hahn spin-echo sequence, (c) CPMG ((j>2 = <j>i) 
and CPMG-2 (<j)2 — cj>i + n) sequences, (d) PDD sequence, 
(e) CDD sequence of order n, CDD„ = C n . (f) UDD sequence 
scheme with 4 pulses, i.e. UDD of order 4, UDD4. 



D. Carr-Purcell (CP) and 
Carr-Purcell-Meiboom-Gill (CPMG) 

To avoid the decay of the echo due to the finite cor- 
relation time of the environment, Carr and Purcell (TT| 
reduced the cycle time by splitting the total time into 
shorter segments of equal length, and a refocusing pulse 
in the middle of each segment. 

Figure HJc) shows the pulse sequence for an initial con- 
dition of S y with <pi = 4>2 = y- The resulting evolution 
operator is f T /2Y f T Y fr/2- Later on, Meiboom and Gill 
|F3 | suggested to shift the phase of the refocusing pulses 
by 7r/2, so that the rotation axis is the same as the orien- 
tation of the initial state. For perfect pulses, both cases 
are equivalent, but only the CPMG version compensates 
flip-angle errors of the refocusing pulses. For a flip-angle 

error exp |— i§ ei 6i >e j = exp |— iSyAuiiTpj for every i 
in Eq. (|23p . the zero-order average Hamiltonian is pro- 
portional to AuiiSy. It thus commutes with an initial 
condition along the y axis, (the CPMG case) and has no 
effect, but it causes an unwanted rotation of an initial 
state oc S x (the CP case). In the following we call this 
sequence with identical 7r-pulses CPMG. 

An alternative sequence that also compensates flip- 
angle errors of the refocusing pulses is shown in Fig. [He), 
with <p2 = <pl + tt = —y- For hard pulses and vanishing 
delays between the pulses, the zero-order average Hamil- 



tonian of this sequence vanishes, for arbitrary flip-angle 
errors, and the first non- vanishing term is of order t c /tb 
and proportional to S x . As a consequence, an initial 
condition proportional to S x is less affected under this 
sequence. In what follows, we will call this DD sequence 
CPMG-2. 

The effect of pulse errors during CPMG and CPMG-2 
on the spin dynamics was studied in Refs. (43 - l50| . and 
we will show some effects in the following sections. 



E. Periodic Dynamical decoupling (PDD) 

A sequence called XY-4 in the NMR community was 
proposed initially to compensate the sensitivity of the 
CPMG-like sequences to non-perfect pulses [ill [p4| . 
Later, it was found equivalent to the shortest universal 
DD sequence that cancels the zero order average Hamilto- 
nian for a general SE interaction of the form [3jl [5l| . 
This sequence, depicted in Fig. [H[d) and called periodic 
dynamical decoupling (PDD), has an evolution operator 
of the form Yf T Xf T Yf T Xf T . Because it suppresses SE 
interactions of the form (|2~4"1) . it compensates errors of 
non-ideal pulses at the end of the cycle. 



F. Concatenated Dynamical Decoupling (CDD) 

The concatenated DD (CDD) scheme pj], HH recur- 
sively concatenates lower order sequences to effectively 
increase the decoupling order. The CDD evolution oper- 
ator for a recursion order of n is given by 



CDD, 



YC n -iXC n -iYC n -iXC n -i, 



(26) 



where C = f T and CDDi = PDD. Fig. HJe) shows a 
general scheme for this process. Each level of concate- 
nation reduces the norm of the first non- vanishing order 
term of the Magnus expansion of the previous level, pro- 
vided that the norm was small enough to begin with. The 
latter reduction is at the expense of an extension of the 
cycle time by a factor of four. 



G. Uhrig dynamical decoupling (UDD) 

Uhrig proposed a different approach to the goal of 
keeping a qubit alive [2 (J . |52j : For a given number N of 
pulses during a total time r c , at what times should these 
pulses be applied to minimize the effect of the system- 
environment interaction? The solution he found for the 
times ti 



is 



h = t c sin 



2(JV+1) 



(27) 



where tjy+i = r c is the cycle time and to — the starting 
time. Defining n = t{ — ti-i the UDD evolution operator 
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for a sequence of N pulses is 

UDDyv = f TN+1 Yf TN Y...Yf T2 Yf T1 



(28) 



and its schematic representation is given in Fig. [^f). 
The CPMG sequence is the simplest UDD sequence of 
order N = 2. 

Cywinski et at explained the performance of the DD 
sequence by finding its spectral filter for the bath-modes 
|5a |. They found that the effect of the UDD pulse se- 
quence leads to an efficient spectral filter for slow mo- 
tion bath-modes. It was shown that UDD is the best 
sequence for reducing the SE interaction in the limit of 
low- frequency noise [20l HH, |53| . Rigorous performance 
bounds for the UDD sequence were found by Uhrig and 
Lidar in Ref. |54|. 



IV. EXPERIMENTAL RESULTS 
A. System and environment 

Experiments were performed on a polycrystalline 
adamantane sample using a home-built solid state NMR 
spectrometer with a 'H resonance frequency of 300 MHz. 
The adamantane molecule contains two nonequivalent 
carbon atoms. Under our conditions, they have simi- 
lar dynamics. Working with natural abundance (1.1 %), 
the interaction between the 13 C-nuclear spins can be ne- 
glected. The main mechanism for decoherence is the in- 
teraction with the proton spins. As discussed in section 
ITT1 this interaction (p~2|) is not static, since the dipole- 
dipole couplings within the proton bath cause flip-flops 
of the protons coupled to the carbon. 

Considering that all the proton spins / are equiva- 
lent, we can estimate the correlation time of the time- 
dependent SE interaction (fT2"j) with the decay time of 
the correlation function i{(t). While the correlation func- 
tions P z (t) of Eq. (TT3")) cannot be measured directly be- 
cause we cannot address individual spins of the bath, we 
get a very good estimate by measuring 



ixit) 



Tr 



Tr 



{t(0)t(*)} 
{4(0)4(0)}' 



(29) 



i.e. the proton free-induction decay (FID) (solid line 
in Fig. [3]). The time evolutions in equation (|29|) 
are determined by the bath Hamiltonian "He, Ix{t) = 



-m E t 



Et I x e %UEt and I x = £). ■ P x . Simulating it with the 
dipole-dipole Hamiltonian of Eq. ([TO and using the same 
Hamiltonian for calculating ii,(t), we find the correla- 
tion function represented with dashed line in Fig. |3] 
The spectral density of the bath is well approximated 
by a normal (Gaussian) distribution and the system- 
environment interaction is weak compared with the intra- 
bath interaction (|6<?j| tb ^$ 1/3). 

The 7r pulses for DD were applied on resonance with 
the 13 C spins. Their radio-frequency (RF) field of 2n x 
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Figure 3: (color online) Evolution of the normalized spin cor- 
relation functions for the bath spins (protons). The solid line 
represents the proton FID signal [i x (t)] and the dashed line 
the numerically simulated i z (t). 



48kHz gives a 7r-pulse length of r p =10.4/is. The measured 
RF field inhomogeneity is about 10%. We performed 
experiments where the delay t between successive DD 
pulses was varied from 10/is to 200/iS. We prepared the 
initial state by using the sequences of Fig. (5Ja) and we 
measured the survival probability of the magnetization 



s u (t) 



Tr 



Tr 



{s u {0)S u {t)} 
{s u (0)S u (0)}' 



(30) 



where u — x,y,z. The solid line of Fig. |4] shows the 
experimental observation of this survival probability from 
an initial condition S x under a free evolution ( 13 C FID). 



B. Hahn echo 

As shown in Fig. U , the decay of the S-spin mag- 
netization is reduced by the Hahn echo sequence. The 
results of the Halm echo are marked by square points. 
Compared to the free induction decay, the decay rate is 
reduced approximately by a factor of 2. 



C. CPMG 

Figure [5] shows the experimental results of the CPMG 
sequence of Fig. [^c) with <pi — (j>2 = V- Different rows 
correspond to different initial conditions S x , S y and S z 
of the 13 C qubit. The left hand panels show the survival 
probability (|30|) as a function of the total evolution time 
t = Mt c (including the pulses), and the right hand panels 
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Figure 4: (color online) Survival probability of the S'-spin un- 
der free evolution ( 13 C FID) and after a Hahn echo sequence. 
An initial condition S x was prepared. 



show the same data as a function of the number of applied 
pulses. 

The plots show that the decay of the survival probabil- 
ity depends crucially on the initial state of the qubit; we 
shall henceforth refer to the initial state in the direction 
of the DD pulses as the "longitudinal" state, and the ones 
perpendicular to the pulses as the "transverse" states. 
Flip-angle errors, which arise from inhomogeneous radio- 
frequency fields, are usually the dominant imperfection 
in this type of experiments. When the CPMG sequence 
is applied to a longitudinal initial condition, flip-angle 
errors do not affect the performance of the decoupling, 
since they are compensated over each cycle consisting 
of two pulses As a result, the decay rates for lon- 

gitudinal states are about an order of magnitude lower 
than for transverse initial conditions. We also observe 
an unexpected oscillation pattern for transverse initial 
states. These kind of strong asymmetries have been re- 
ported in different samples, and have been hypothesized 
to be due to stimulated echoes induced by pulse errors 
[3 HE EE or due to the non- negligible effects of the 
interaction Hamiltonian acting during the finite width 

pulses IHEzim. 

The right hand panels show the same data, but plotted 
against the number of pulses. They clearly show that the 
oscillation frequency depends on the number of applied 
pulses or equivalently on the total pulse-irradiation time. 
Similar oscillations have been also reported in different 
samples (45l [47- 49] . In our experiments, the oscillation 
pattern originates from the bimodal distribution of rf 
field amplitudes in the coil. For our present analysis, the 
beating is not important because it could be reversed [48j 
or avoided by improving the rf field coil. We instead con- 
centrate on the decay of the envelope, which represents 
the overall survival probability of the signal. 
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Figure 5: (color online) Magnetization evolution for the 
CPMG sequence. From top to bottom the initial conditions 
are S y , S x and S z . The left-hand panels represent the C 
magnetization as a function of the total evolution time while 
right-hand panels show its evolution as a function of the num- 
ber of applied pulses. The legend at the bottom gives the 
delays r between successive pulses. 



For the longitudinal initial state (upper panel), panel b 
shows that the signal decay, as a function of the distance 
r between successive pulses, remains constant until r — 
30/iS. The corresponding cycle time is t c = 2t + 2t p = 
80.8/xs, which is comparable to the bath correlation time 
tb\ hence; the signal decay for cases below r = 30/is is 
mainly due to pulse errors. For longer delays r c > tb 
(jb ~ 110/is), the decay rate increases because of the 
reduction of time reversal efficiency in the fluctuating 
environment. 



D. CPMG-2 

Figure [6] shows the corresponding results for the 
CPMG-2 sequence. Since the first non-vanishing order 
of the Magnus expansion for the CPMG-2 sequence (con- 
sidering flip-angle errors) commutes with S Xl we expect 
that the signal decay for the S x initial state is similar 
to that of the longitudinal initial state of the CPMG ex- 
periments. The experimental results shown in (Fig. |6]) 
clearly agree with this expectation. For the other ini- 
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Figure 6: (color online) Magnetization evolution for the 
CPMG-2 sequence. From top to bottom the initial condi- 
tions are S x ,S y and S z . The left-hand panels represent the 
13 C magnetization as a function of the total evolution time, 
while right panels show its evolution as a function of the num- 
ber of applied pulses. The legend at the bottom gives the 
delays r between successive pulses. 



vs. y. The decays still contain a small oscillatory contri- 
bution. Since it appears to depend mostly on the num- 
ber of pulses, rather than on the delays between them, 
we attribute them to pulse errors that are not completely 
canceled. Compared to CPMG and CPMG-2, the period 
of the oscillation is one order of magnitude longer, indi- 
cating that the effect of the pulse imperfections has been 
reduced by an order of magnitude. This general improve- 
ment against pulse errors is because the sequence cancels 
the zeroth order average Hamiltonian of the more general 
SE interaction (|24|) . while the CPMG sequences cancel 
only its pure-dephasing part. 

In Fig. [T](b) and (d), the decay rates, in units of pulses, 
up to r = AOfis, i.e. t c = 4 (r + t p ) — 201.6^s, are 
equal to within experimental error. This shows that the 
sequence is more robust, compared to previous sequences, 
in the regime where r c exceeds the bath-correlation time. 
However, from the time evolution of the left panels a and 
c, the decay rate is larger than the longitudinal case of 
CPMG or the S x case of CPMG-2. 

If the initial state is proportional to S z , its evolution is 
qualitatively different. We believe that this results from 
the fact that it is parallel to the static field and com- 
mutes with the free precession Hamiltonian ([5]). As a 
consequence, this evolution reflects the implementation 
errors of the sequence. The source of the decay is the 
pulse errors due to which the average Hamiltonian no 
longer commutes with the initial state. This experiment 
provides a means of quantifying pulse errors, and cali- 
brating an optimal setup of the sequence to enhance its 
performance. 



CDD 



tial states, an oscillatory behavior similar to that for the 
transverse state of the CPMG is observed. Although the 
oscillation still depends on the number of pulses (right 
panels), the frequency is slower than in the CPMG case, 
and the envelope of the oscillations decays more slowly. 
The origin of the oscillation pattern is again the bimodal 
distribution of the inhomogeneity of the RF field generat- 
ing an effective field along the x axis. The experimentally 
observed oscillation agrees with the results of the effec- 
tive nutation experiment. Again, we will concentrate on 
the decay of the envelope. 



The qualitative behavior of the CDD experiments is 
similar between different orders and to the PDD one, but 
they change in the time scale for which the initial state 
can be maintained. They are also more robust against 
pulse errors - the oscillation pattern is not observed. For 
details of the survival probability evolution see the ap- 
pendix [A] A summary of the results is presented in Fig. 
|8]in section IVl where the decay times for different CDD 
orders and delays r between pulses are plotted. 

G. UDD 



E. PDD 

Figure [7] shows the signal decay for different initial con- 
ditions under the application of the PDD sequence of Fig. 
[2^d). One observes that the signal decay evolves quali- 
tatively similar for initial conditions in the plane trans- 
verse to the static field, i.e. S x and S y . This agrees with 
the theoretical predictions: the sequence of evolutions 
Y f T X f T Y f T X f T is nearly symmetric with respect to x 



The experimental survival probabilities for the UDD 
sequences have the same qualitative behavior as the 
CPMG curves. They manifest the same asymmetries 
with respect to the initial state. That is expected because 
UDD sequences also only reduce SE interactions of the 
form (fTTTj) . We observed that with increasing UDD order, 
the decay rates increase, as predicted in Ref. [54j for the 
conditions satisfied in our experiments where \bsj Tg < I 
and t c ~ tb . A summary of the rates is shown in the 
next section in Fig. [5J An extensive analysis of the per- 
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Figure 7: (color online) Signal decay of the initial state of the 
qubit for the PDD sequence. From top to bottom the initial 
conditions are S x ,S y and S z . The left-hand panels represent 
the decay as a function of the total evolution time while the 
right-hand panels show the decay as a function of the number 
of applied pulses. The legend at the bottom gives the delays 
r between successive pulses. 



formance of UDD sequences and non-equidistant pulse 
sequences against equidistant ones for the present exper- 
imental conditions will be given elsewhere (55j . 



V. COMPARISONS: OPTIMAL CHOICES 

The goal of DD is the preservation of quantum states 
by the application of suitable decoupling sequences. If 
the pulses are ideal and they are applied with very 
short delays, it is possible to preserve quantum states 
for arbitrarily long times in the presence of a system- 
environment coupling that is linear in the system oper- 
ators. However, for experiments using non-ideal pulses, 
a finite cycle time optimizes the DD performance. This 
can be seen very clearly in the summary of the experi- 
mental results presented in Fig. [5] The left panels show 
the DD decay times as a function of the delay r for every 
sequence and different initial conditions. As an example, 
for CPMG when the initial condition is longitudinal to 
the pulses (S y ), the optimal cycle time is t c = 80.8/is 



Figure 8: (color online) Relaxation times of different initial 
conditions under DD conditions as a function of the delay 
between pulses r (left panels) and the cycle time r c (right 
panels). From top to bottom the initial condition is given by 
Po = S x , S y and S z respectively. An optimal r and conse- 
quently t c is observed for each sequence. The reduction of 
the relaxation time to the right side of the optimal value is 
due to the shifting environment: in this regime the cycle time 
is longer than the correlation time of the bath, r c > tb- The 
reduction for short cycle times indicates that in this regime, 
accumulated pulse errors dominate. 



(t = 30/j.s). For longer cycle times, the decay time gets 
shorter, since the environment changes during the cycle 
and the refocusing efficiency decreases. 

While shorter cycle times should give even better re- 
sults under ideal conditions, we find experimentally a de- 
crease of the relaxation time. This can be attributed to 
an accumulation of pulse errors, which dominates in this 
regime. Similar results are observed for the CPMG-2 if 
we exchange S y with S x . This can be seen clearly un- 
der conditions of transverse initial states (S x for CPMG) 
where the decay time is proportional to the cycle time. 
This means that the error per cycle is independent of 
the cycle time and corresponds thus to a zero-order term 
of the average Hamiltonian. This is the behavior ex- 
pected for flip angle errors, which are the main source of 
the decay in this regime. Since flip angle errors are in 
no way compensated for transverse states in the CPMG 
sequence, their accumulated effect is so strong that the 
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Figure 9: (color online) CDD order n as a function of its 
optimal delay between pulses T opt to reduce decoherence. The 
experimental square points seem to satisfy a relation given by 
n — c — bln(r op t/r_e). The solid (red) line shows a fitting 
curve with parameters c = (0.9 ± 0.2) and b = (—0.9 ± 0.1). 



optimal cycle time exceeds tb and the sequence performs 
only marginally better than the Hahn echo sequence, 
which has the longest cycle time. 

If we consider the CPMG-2 sequence with the initial 
conditions S y and S z , the decay time grows oc t\ for short 
times. This implies that in this case, the dominant error 
term is proportional to r c , i.e. it corresponds to a first- 
order term of the average Hamiltonian. Moreover, its 
optimal relaxation time is one order of magnitude longer 
than the Hahn echo decay time. 

The behavior of the UDD sequences is similar to that 
of CPMG. We show here only their decay times for S y as 
initial condition [Fig. [SJc)], i.e. longitudinal to the DD 
pulses. They are plotted as a function of the average de- 
lay between pulses. The figure shows that the UDD decay 
times are always shorter than those of CPMG. Moreover, 
increasing the UDD order reduces the decay time, as ex- 
pected by theoretical expectations when \bsj \ tb < 1 and 
t c ~ Tb 0| . A regime where UDD performs better than 
CPMG may perhaps exist at short cycle times compared 
with tb [55J, provided the pulse errors can be made suf- 
ficiently small that they do not dominate over external 
sources of decoherence. Recent proposals of UDD based 
sequences that reduce decoherence of a general SE inter- 
action like Eq. (|24f . could allow one to find this regime 

nam. 

For PDD the optimal cycle time for S x and S y is 
t c = 321.6/is, (r = 70/is), which is longer than rg. The 
resulting performance is relatively poor. The resulting 
decay time is similar to that of CPMG for the same cycle 
time [Fig. ltd)] or CPMG-2 for the S x initial condition 
[Fig. H{b)]- However, for these particular initial condi- 
tions CPMG and CPMG-2 can be made to perform an 
order of magnitude better by reducing the cycle time. 



For CDD2 with initial conditions transverse to the 
static field, the optimal cycle time is r c = 16r + 20r p = 
688/is, i.e. r = 30/zs. For CDD 3 and CDD 4 the shortest 
delay time between pulses of t = 10/is and r = 2.5/^s are 
the optimal situations, giving t c (CDD3) = 64t + 84r p = 
1513.6/xs and t c (CDD 4 ) = 256t + 388t p = 4675. 2^s. The 
optimal delay t becomes shorter with increasing CDD or- 
der, because the cycle time increases by a factor of 4 for 
each level of concatenation. Apparently, the pulse errors 
do not accumulate as strongly as in the case of CPMG, 
which may be attributed to the fact that CDD is designed 
to compensate pulse errors (l7l . |33T | . The crossover cycle 
time, where the transition occurs from a decay dominated 
by pulse errors to the regime where the decay is domi- 
nated by the short bath correlation time is increased, as 
shown in the right hand panels. Even for cycle times that 
are much longer than the bath-correlation time, the CDD 
provides a significant reduction of the decoherence rate 
compared to the free evolution decay (dotted lines) and 
the Hahn echo decay (dashed lines). An optimal cycle 
time when considering imperfect finite pulses was pre- 
dicted by Khodjastch and Lidar [33]. Figure [5] shows the 
experimental relation between the optimal delays r and 
their respective CDD order n (square points). It seems 
to satisfy a relation given by n = c — b\n(r op t / tb) , where 
c and b are constants and r op t is the optimal delay for a 
given n (see below). 

The unifying result of the curves shown in the left- 
hand panels of Figure|8]is that the optimal delay between 
pulses is always shorter than the bath correlation time, 
with comparable values for all sequences, with the single 
exception of the CPMG sequence for initial conditions S x 
and S z , as discussed above. Clearly, this timescale is de- 
termined by the (average) delay between pulses r, not by 
the cycle time r c . Expressing this differently, one might 
say that only a small fraction of what is lost in a single 
echo can be refocused by compensated sequences. If we 
look at pulse spacings longer than the bath correlation 
time tb, the differences between sequences become very 
small and the decay times approximate those of FID and 
Hahn echo. Accordingly, it appears important to keep 
the number of pulses per unit time constant when com- 
paring different DD sequences. 

We now compare the different DD sequences with the 
optimal cycle time for each sequence. Figure [TU] shows 
the evolution of the survival probabilities for different 
initial conditions for all the sequences discussed here. 

As a general rule, we note that for increasing CDD or- 
der, the optimal pulse delay r gets shorter and for longer 
delays between pulses, higher CDD orders do not perform 
better than lower CDD orders. Hence, keeping the delay 
between pulses constant, there is an optimal CDD order 
for reducing decoherence as predicted in Refs. [33j,|34|. It 
is difficult to find accurately the optimal CDD order as a 
function of r from Figs. [8] (a) and (c) in order to com par e 
with the theoretical predictions of Eq. (140) in Ref. [3J|. 
However, the relation given in Fig. [5] for the optimal de- 
lay between pulses T opt behaves similar. A linear fitting 



12 



of the experimental data gives n = c — 61n(r op t / Vb) with 
c = (0.9 ± 0.2) and b = (-0.9 ± 0.1) agreeing well with 
the predicted expression [34J ] . 

The best DD sequence and its corresponding optimal 
cycle time depends on hardware limitations, and impor- 
tantly on the desired goal. If one aims to freeze a quan- 
tum state during a short time, its value will bound the 
cycle time and as a result, the maximal CDD order that 
can be applied. For longer times, increasing the CDD 
order will be advantageous, but power dissipation may 
force a reduction in the number of pulses and simulta- 
neously the CDD order. For specific initial conditions, 
CPMG and CPMG-2 are the best choices for reducing 
decoherence; however the large asymmetry of these se- 
quences to other initial conditions limits their usefulness 
when the initial state of the qubit is unknown. Note that 
in these cases, it has been shown that coherences could 
be frozen as labelled polarization [50(. If the goal is the 
preservation of an unknown quantum state, the CDD se- 
quences provide the best overall performance. 

While the SE interaction produces pure dephasing to 
the spin, in principle DD sequences that compensate pure 
dephasing decoherence should be sufficient. In conse- 
quence concatenating sequences like f T Yf T Y, which re- 
duce pure dephasing processes, could be beneficial. How- 
ever, the finite precision of control pulses generates an ef- 
fective Hamiltonian of the form (f!H)) . and thus sequences 



developed to reduce pure dephasing processes have asym- 
metric performances against initial states directions, i.e. 
they do not generate a unit evolution operator of the 
qubit. Thus, we compared, as a test bed, CPMG/UDD 
sequence that compensate pure dephasing with XY-4 
[l3L [Tij based sequences that compensate a general inter- 
action in order to show their effects against pulse errors. 

From our results it is evident that for short delays 
the main source of DD decays are static pulse errors; 
this limits the maximal performance. However, CPMG 
and CPMG-2 show the potentially achievable DD perfor- 
mance if the pulse errors are reduced. This implies that 
new DD proposals should focus on the compensation of 
pulse errors for these kinds of experimental conditions, 
similar to the prop osal by Viola and Knill [l5j or Uhrig 
and Pasini [56^ [H3|- CDD- type sequences do compensate 
for pulse errors, but only at the end of the CDD cycle. 
This limits their performance because of the exponen- 
tial growth of the cycle time with the CDD order. As 
an alternative method, we suggest to find sequences that 
compensate pulse errors to zero order during each step 
of the concatenation procedure. That would be advan- 
tageous because the zero order compensation cycle time 
remains constant and equal to the PDD cycle time, as 
assumed for ideal pulses. 



VI. CONCLUSIONS 
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Figure 10: (color online) Time evolution of the survival prob- 
ability for the optimal cycle times of the different DD se- 
quences. From top to bottom the initial states are S x ,S y and 
S z . The optimal delays r are given in the legends. 



We have experimentally applied different dynamical 
decoupling sequences to a qubit-system coupled to a spin- 
bath in order to test and compare their performance. The 
system used is typical for spin-based solid-state systems 
where the spectral density of the bath is given by a nor- 
mal (Gaussian) distribution and the system-environment 
interaction is weak compared with the intra-bath inter- 
action. The experiments were performed in the regime 
where the average spacing between the pulses is compa- 
rable to the bath-correlation time. This article focuses on 
measuring and fighting decoherence, and the results do 
not depend on the readout or initialization scheme used 
for that purpose. Thus, the results should apply directly 
to other spin-based quantum information processing sys- 
tems, such as electron spins in diamonds [30h32| |. electron 
spins in quantum dots [H, [2^, |38j and donors in silicon 

While the design of DD sequences is typically based 
on the assumption that the cycle times is shorter than 
the bath-correlation time, we demonstrated that even 
without satisfying this condition, dynamical decoupling 
reduces decoherence significantly. We showed that the 
main limitation to the reduction of the DD decay rates 
is due to the finite precision of the control operations 
- in our system, flip-angle errors were the main source. 
Therefore, CPMG or UDD-type sequences that reduce 
purely-dephasing or purely-dissipative interactions with 
the bath perform well only for specific initial conditions 
of the qubit ensemble. For the privileged initial condi- 
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tion, CPMG-type sequences performed better than any 
other DD sequence. But, if the goal is to approach a 
unit evolution operator, PDD sequence and its concate- 
nated form (CDD) are the best overall option. In agree- 
ment with previous predictions [33j-|35j our results show 
that, depending on limitations of hardware and the re- 
quired time to keep the initial state coherent, increasing 
the CDD order is not always useful. There is an optimal 
CDD order depending on the power available for the con- 
trol pulses and their finite precision. We present strong 
evidence that in order to improve dynamical decoupling 
sequences, they should be designed to compensate pulse 
errors. 
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Figure 13: (color online) Signal decay of the initial state of the 
qubit for the CDD4 sequence. From top to bottom the initial 
conditions are S x ,S y and S z . The left-hand panels represent 
the decay as a function of the total evolution time while the 
right-hand panels show the decay as a function of the number 
of applied pulses. The legend at the bottom gives the delays 
r between successive pulses. 



